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Abstract: The challenge of demonstrating that the matter produced in heavy ion colli- 
sions is a deconfined quark-gluon plasma, as predicted by lattice QCD calculations, is the 
challenge of measuring the number of thermodynamic degrees of freedom v ~ e/T 4 at the 
time to at which the matter comes into approximate local thermal equilibrium and begins 
to behave like a hydrodynamic fluid. Data from experiments done at the Relativistic Heavy 
Ion Collider have been used to estimate to and to put a lower bound on the energy density 
e(to)- However, measuring v has seemed out of reach, because no current data serve even 
as qualitative proxies for the temperature T(io)- We point out that v may equally appro- 
priately be defined via v ~ s 4 /e 3 , where s is the entropy density, which can be estimated 
from the measured final state entropy. This estimate is based on the testable assumption 
of an isentropic expansion. The observation of jet quenching has the potential to provide 
an upper bound on the energy density at early times. Our goal is to motivate such an 
analysis by pointing out that it would set a lower bound on v. 
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Measurements of the final state produced in heavy ion collisions at */sj^ = 1^0 and 
200 GeV at the Relativistic Heavy Ion Collider (RHIC) in Brookhaven have provided strong 
evidence that the matter created in these reactions is rapidly thermalized and has many 
properties expected from a strongly coupled quark-gluon plasma What is lacking, 
however, is direct evidence that could be linked to the prevalence of unconfined colored 
excitations, which are thought to be the defining feature of this novel state of matter. 

At present, the best evidence for deconfinement of quarks may be that derived from the 
phenomenological success of valence quark recombination models [|J, which are based on 
the assumption of the existence of a thermalized phase of unconfined "constituent quarks" 
immediately prior to hadronization. In particular, the quark recombination models explain 
the different saturation levels and thresholds for the elliptic flow of different hadrons in 
terms of a universal elliptic flow pattern at the valence quark level. A serious limitation 
of this evidence is that it addresses physics near the quark-hadron transition, which is 
inherently nonperturbative and thus not amenable to controlled theoretical approximations. 
It would be desirable to obtain evidence for the liberation of colored degrees of freedom in 
the matter created at RHIC which relates to thermalized matter present during an earlier 
stage of the reaction, when the temperature is higher, and to observables that allow for a 
more controlled theoretical description. 

In QCD with quark masses as given in nature, there is no difference in symmetry 
between an equilibrated gas of hadrons and a quark-gluon plasma, implying that the tran- 
sition between these two regimes of strongly interacting matter may be a continuous, albeit 
possibly rapid, crossover ||. Indeed, this is what is predicted by ab initio calculations of 
QCD thermodynamics done using the methods of lattice gauge theory Hence, an ex- 
perimental demonstration of deconfinement cannot be seen as the answer to some "yes/no" 
question, but must instead involve the measurement of some physical property of the matter 
created in heavy ion collisions that can also be predicted by controlled theoretical calcula- 
tions, and which takes on quite different values below and above the crossover between a 
hadron gas and a quark-gluon plasma. 

Because QCD is an asymptotically free theory, there are many quantities whose cal- 
culation is controlled at temperatures far above the crossover, where the effective QCD 
coupling a s (T) becomes small. However, we do not expect that collisions at RHIC can 
create matter in this asymptotic regime. This means that in looking for quantities whose 
theoretical calculation is under control, we must ask what quantities can be calculated by 
rigorous numerical methods. Lattice calculations of non-thermodynamic quantities, like 
spectral functions || and viscosities §, are still in their infancy and are presently re- 
stricted to QCD in the quenched approximation. In contrast, lattice calculations of QCD 
thermodynamics have reached maturity because they can be formulated conveniently in 
Euclidean quantum field theory, the natural arena for lattice QCD. Lattice simulations of 
the QCD equation of state including dynamical quarks of various numbers of flavors are 
available, which permit model- independent conclusions Q. 

In QCD thermodynamics, there is one observable whose value changes by more than an 
order of magnitude across the crossover transition from hadron gas to quark-gluon plasma, 
namely the effective number of thermodynamic degrees of freedom v(T). A common defi- 
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nition of v{T) is via the relation 



e{T) = ^v{T)T\ (1) 



where e and T are the energy density and temperature, respectively. For an ideal gas of 
massless, noninteracting constituents, v counts the number of bosonic degrees of freedom 
plus the number of fermionic degrees of freedom weig hted by 7/8. Eq. @ makes the mea- 
surement of v in heavy ion collisions seem a remote possibility, because there is nothing in 
the current suite of data from RHIC that is thought to serve as a proxy for the temperature 
at early times. The temperature at freezeout, when the matter is again hadronic, is well 
determined from many measurements of the hadronic final state. A measurement of the 
temperatures of the quark-gluon plasma presumed to be present at early times is one of the 
goals of studies of direct photon and dilepton emission in heavy ion collisions but, so far, 
no evidence for thermal photon or dilepton radiation has been observed at RHIC. Even if 
they will be observed eventually, thermal photons will only yield a weighted time-average 
of the temperature, which may not be sufficient to determine the function v{T). 

Here, we make the simple observation that v can equally well be defined via the entropy 
density 

2tt 2 

s(T) = —u(T)T 3 , (2) 

and hence via 

For an ideal gas of massless degrees of freedom, the expressions (||), @ and (||) are equiv- 
alent definitions of v. We shall take (||) as our definition, because in so doing we realize 
that v{T) can be measured without measuring T itself. An ideal gas of ultrarelativistic 
pions has v = 3, whereas an ideal gas of noninteracting gluons and three flavors of massless 
quarks has v = 47.5. 

Lattice QCD calculations show that v increases by more than a factor of ten over a 
narrow range of temperatures centered at a crossover temperature T c = 170 ± 10 MeV Q. 
This rapid increase in the number of degrees of freedom is the direct consequence of de- 
confinement at high temperatures in QCD and, if we limit ourselves to thermodynamic 
observables, the measurement of v is the only possible route to an experimental demon- 
stration of deconfinement. 

QCD does not describe an ideal gas of noninteracting quarks and gluons, except at 
infinite temperatures. And indeed, there is growing evidence from the RHIC experiments 
that the matter they are creating is strongly interacting, thermalizing rapidly, flowing like 
a liquid, and opaque to energetic partons produced within it In the case of strongly 
interacting matter, the different definitions of v provided by ([[]), (Q) and (Q) are not equiva- 
lent, and furthermore v cannot be considered as a count of well-defined degrees of freedom. 
Indeed, in an ideal liquid there are no well-defined long-lived quasiparticles. Neverthe- 
less, upon choosing one definition — and we shall choose (||) — v is reliably calculable 
on the lattice, potentially measurable as we shall discuss, and remains a valid measure of 
deconfinement . 
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Figure 1: The effective number of degrees of freedom of a hadronic resonance gas, using the three 
definitions of v provided by Eqs. (1-3): (1) long-dashed line (blue), (2) short-dashed line (green), 
and (3) solid line (red). The solid curve moves by about 5% depending on whether we use all 
established resonances in the particle data book or only those commonly included in the chemical 
freezeout analysis. This can be considered an estimate of the theoretical uncertainty in the curve. 

Lattice QCD calculations indicate that throughout the temperature range 2T C < T < 
5T C , v(T) is between 70% and 80% of that for an ideal quark-gluon plasma ||], meaning that 
33 < v < 38. At lower temperatures, closer to the crossover, the value of v extracted from 
lattice QCD calculations is more significantly different for different definitions of v, because 
just above the crossover the pressure deviates more from its ideal gas value than the energy 
density does. At T = 1.5 T c , for example, with the definition ([[]) lattice calculations § 
yield v ~ 37, whereas with our definition @ they yield v ~ 27. 

There is an additional benefit to defining v via (|3|) as we do, over and above making 
the quantity observable in practice. In order to be convincing, a future experimental 
measurement of v that agreed with what lattice QCD predicts for the quark-gluon plasma 
would need to discriminate between this prediction and that of a hadronic resonance gas 
with T < T c . We see from Fig. 1 that, with our definition of v, this requires excluding 
v ~ 10. The three different definitions of v all agree for a gas of noninteracting massless 
particles, but massive particles — as in a resonance gas — contribute somewhat less to 
s/T 3 than they do to e/T 4 , and hence less still to s 4 /e 3 . This makes our definition of v 
advantageous. 

Before turning to the challenges associated with measuring s and e, we note that the 
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value of v defined from them via (||) can only be compared to QCD thermodynamics if 
there is independent evidence that the matter under study is in (approximate) local thermal 
equilibrium. At RHIC, such evidence is believed to be provided by the agreement of the 
elliptic flow measured in noncentral collisions with hydrodynamic model predictions Q. 
Such predictions are based on the assumption that the matter behaves like a fluid in local 
thermal equilibrium, with arbitrarily short mean free paths and correspondingly strong 
interactions. The fact that the data show as much elliptic flow as predicted indicates 
that this assumption must already be valid soon after the collision, early enough that the 
azimuthal spatial anisotropy due to the nonzero impact parameter has not had time to be 
significantly reduced via free streaming of weakly interacting quasiparticles. Quantitatively, 
it is estimated that local thermal equilibrium and the onset of hydrodynamic behavior must 
occur by a time to ~ 0-6 — 1-0 fm/c Q. 1 We shall take to = 1.0 fm/c in the following. In 
the remainder of this paper, we shall discuss how s(io) and e(io) may be estimated, using 
present and near-future data. 

The estimate of s(io) relies on the fact that the entropy of an ideal fluid is conserved 
during its hydro dynamical evolution. The value of s(to) can therefore be deduced from 
the value of dS/dy at chemical freezeout, upon assuming that the fluid remains in local 
thermal equilibrium between the time to when this condition is established and the time 
at which chemical freezeout occurs. Note that we shall need to know the volume of the 
system at time to, which is early enough that little transverse expansion has taken place 
making this easy to estimate, but we shall not need to know the volume of the system at 
freezeout. Indeed, we do not need to know anything about the system at freezeout except 
its entropy. In an isentropic expansion, it is entropy that is conserved, regardless of how 
entropy density evolves. 

There are, in principle, two ways of estimating dS/dy. One approach uses an analysis 
of the composition of the fireball at chemical freezeout to derive the entropy per hadron, 
which can then be used to relate the measured charge particle multiplicity dN c ^/dy to the 
entropy dS/dy, assuming full thermal phase space occupation at the freezeout time. The 
other approach Q uses the measured multiplicities of stable hadrons together with experi- 
mental data on Hanbury-Brown-Twiss (HBT) two-particle interferometry to estimate the 
full phase space distributions fi{r,p) at kinetic freezeout, and from them the entropy 

S = T.j ^^[-/<ln/ 4 ±(l±/,)ln(l±/ 4 )] , (4) 

i 

where the upper (lower) sign holds for bosons (fermions). 

We begin with a review of the results of the kinetic freezeout analysis. Pal and Pratt || 
have used final state multiplicities and "radius" parameters deduced from a HBT corre- 
lation analysis to determine dS/dy in data from = 130AGeV Au+Au collisions. For 
collisions in a 11% centrality window, Pal and Pratt find dS/dy ~ 4450 with an estimated 
systematic error of ±400. Two effects may affect the reliability of this estimate. First, the 

1 Recent work suggests that this rapid approach to local thermal equilibrium may occur via plasma 
instabilities, not collisions We caution that if thermal but not chemical equilibrium has been achieved 
at to, then v(to) may be lower than the value calculated in lattice QCD. 
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conventional extraction of HBT radius parameters is based on the assumption that hadrons 
originate from a source that is Gaussian in both position and momentum space. This is 
not a good assumption for a rapidly expanding source, and the measured HBT radius pa- 
rameters do not actually serve as estimates of the radii of the fireball at freezeout. Hence, a 
more direct measurement of the phase space distributions themselves, without introducing 
a parametrization in terms of "radii" , would be of value. And, second, the entropy of the 
hadronic gas must increase as it cools after chemical freezeout but before kinetic freezeout. 
As the mean free path becomes larger and eventually approaches the size of the system at 
kinetic freezeout, the viscosity of the hadron gas grows, ideal hydrodynamics ceases to be a 
good approximation, and entropy is produced. Furthermore, after chemical freezeout more 
resonances decay than are produced, again increasing the entropy. Hence, the entropy at 
chemical freezeout, which is what is of interest to us, must be smaller than that estimated 
at kinetic freezeout. 

Next, we consider the analysis based on the inferred chemical freezeout multiplicities 
of hadrons, which has the potential to extract the entropy at chemical freezeout directly. 
The ratios of hadron multiplicities in RHIC collisions with y/s = 200 AGeV are quite well 
described by the assumption that chemical freezeout occurs at T c h = 170 ± 10 MeV from 
an equilibrated hadron gas at this temperature [1C]. An ideal gas of all established meson 
and baryon resonances at T = 170 MeV has an entropy of about 7.25 per hadron 11, 12], 
which can be compared to the value of 3.6 for an ideal gas of massless pions. After all the 
resonances decay the multiplicity of charged hadrons in the final state is 1.04 per hadron in 
the equilibrated hadron gas 111]]. In the 6% most central collisions with y/s = 200 AGeV, 
the multiplicity of charged particles in the final state is dN^/drj = 665 ± 26 at mid- 
rapidity fi"3[l , corresponding to a dN^/dy that is about 10% larger Jlifl . Putting the pieces 
together, we estimate dS/dy = (665 x 1.1 x 7.25)/1.04 5100 at chemical freezeout at mid- 
rapidity in RHIC collisions with y/s = 200 AGeV. There is a 4% uncertainty in this estimate 
coming from that in the experimental measurement of dN^/drj. The largest theoretical 
uncertainty is that in the factor S/N = 7.25. Changing the chemical freezeout temperature 
by ±10 MeV changes S/N by ±3%. To get a sense of the other uncertainties in S/N, we 
recalculated it using abundances obtained including the widths of states, and found that 
S/N increased to 7.58. There are still further sources of theoretical uncertainty that are 
harder to estimate like, for example, that due to our neglect of resonances not found in the 
particle data book. We do not think this is a large effect, because reducing the number of 
resonances that we include does not have a large effect. We estimate S/N = 7.25 ± 6%, 
and expect that a more systematic analysis could reduce this uncertainty by a factor of 
two. Adding the theoretical and experimental uncertainties in quadrature and rounding 
upwards yields dS/dy = 5100 ± 400. 

We can also compare our estimate of dS/dy at chemical freezeout to that obtained 
by Pal and Pratt at kinetic freezeout. Applying our argument to the 11% most central 
collisions with y/s = 130 AGeV, for which dN^/dr] = 526 ± 20 |l~3| l, yields the estimate 
dS/dy = 4035 ±300. As expected, the central value of this estimate is below that obtained 
by Pal and Pratt, but the difference is within error bars. This indicates that the entropy 
release between chemical and kinetic freezeout is not dramatic. 
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To summarize, the best estimates available at present suggest that dS/dy « 5100 ±400 
at chemical freezeout at mid-rapidity in central RHIC collisions with yfs = 200 AGeV. It 
seems to us that a more careful theoretical analysis can reduce the theoretical uncertainty 
by a factor of two, which would yield a measurement of dS/dy with 5% errors. 

To the extent that the expansion between to = 1 fm/c and chemical freezeout is 
isentropic, the argument first proposed by Bjorken can be used to turn dS/dy into 
a lower bound on s(io). The PHOBOS version of this argument [jlij can be phrased as 
follows. The charged particle multiplicity is reasonably independent of the rapidity y within 
the range \y\ < 1, and all this entropy (S = IdSjdy « 10200 ± 800) must have come from 
within a volume of size 2to7rR 2 at time to- Taking R = 7 fm, this yields s(to) > 33±3 fm -3 . 
The entropy density at to is somewhat greater than this, because of the contribution from 
particles outside \y\ < 1. 

In order to obtain a lower bound on z^(to) using (|3|), it is important to have a lower 
bound on s(to), as the Bjorken argument provides upon assuming isentropic expansion. It 
is crucial, however, to test this assumption. That is, it is crucial to rule out a significant 
increase of the entropy between equilibration at time to and chemical freezeout. Entropy 
production during the hydrodynamic expansion requires some specific mechanism such as 
a strong first order phase transition which can drive the matter away from local thermal 
equilibrium. Lattice QCD calculations indicate that the transition is a crossover, with 
hadronization occuring continuously, but it would be desirable to have experimental con- 
firmation of the absence of a strong first order phase transition, in order to complete the 
justification of our use of dS/dy at chemical freezeout to obtain a lower bound on s(to). A 
strong first order phase transition would lead to large (and possibly non-Gaussian) event- 



by-event fluctuations at low pj- [16]. There is no evidence for such fluctuations in current 
data [O], but given the importance of this issue a more stringent investigation of low-px 
event-by-event fluctuations is called for, looking at the fluctuations of several observables 
and focusing on px significantly smaller than the mean. 

Just as it is important to have a lower bound on s(to), if we wish to obtain a lower 
bound on v we need an upper bound on e(to). The analogous "Bjorken argument", applied 
to dEx/dy (where Ex is the total transverse energy of the hadrons in the final state) only 
yields a lower bound on e(to) because the longitudinal expansion subsequent to to reduces 
dEx/dy. For this reason, we cannot use dEx/dy for our purposes. Putting this another 
way, Bjorken arguments applied to dN/dy and dEx/dy yield lower bounds on both s(to) 
and e(to), and so by themselves these arguments give no constraint on v. 

The analysis of jet quenching data has the potential to yield an upper bound on the 
energy density at early times, as we now discuss. High energy partons traversing strongly 
interacting matter lose energy mainly by gluon radiation after interactions with colored 
constituents of the medium. The theory of this mechanism is well developed within the 



framework of perturbative QCD [Qq, [L9L EG, 21, 22, M. In the multiple soft scattering 



limit, the effect of the medium is encoded in the parameter 

q = p f q 2 dq 2 (da/dq 2 ) , (5) 
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where p is the density of scattering centers in the medium, q 2 denotes the momentum 
transfer in scattering, and da/dq 2 is the differential cross section for scattering of the hard 
parton on a single center. Coherence effects suppress the emission of gluons with energies 
above u c = qL 2 /2, where L denotes the length of material traversed. The mean energy 
lost by the parton is given by AE (3a s CR/7r)u; c , where Cr is the Casimir operator 
for the color representation of the parton. Similar results are obtained in the opposite 
(low opacity) limit, when the interaction with the medium is dominated by a single or few 
scatterings p2fl . In these analyses, the properties of the matter being probed arise only in 
the transport coefficient q. In this sense, jet quenching can be thought of as "measuring 

The relationship between q and the energy density e has not been determined in gen- 
eral, but it is known for the limiting cases of cold nuclear matter [20|, an ideal pion gas, 
and an ideal weakly-interacting quark-gluon plasma . When a high energy parton pen- 
etrates strongly interacting matter, it resolves the partonic constituents of the medium. 
In the case of cold nuclear matter, it mainly interacts with the gluon component of the 
nucleons, and the parameter q can be expressed in terms of the gluon distribution in the 
nucleon [20]. For a medium such as a weakly interacting quark-gluon plasma, in which 
all partons are deconfined, each quark and gluon contributes to the density of scatterers 
independently, and q is directly given by Eq. (||) in terms of the gluon density and per- 
turbatively screened parton-parton cross section. The resulting energy loss in a thermal 
quark-gluon plasma has been calculated by Baier as function of the energy density of the 
plasma [E4j. Remarkably, when the energy loss coefficient is calculated for a thermal gas 
of pions, one finds the same value of q as for a thermal quark-gluon plasma with the same 
energy density. With the benefit of hindsight, this is not entirely surprising, although the 
precision of the agreement may be coincidental. About half of the momentum of a fast 
moving hadron is carried by gluons, and roughly half of the energy density of a weakly 
interacting quark-gluon plasma is contained in gluons. Normalized to the energy density 
of the medium, a fast moving parton can therefore be expected to encounter roughly the 
same number of gluons on which it can scatter. 

It is presumably naive to think that a measurement of q is a measurement of the energy 
density e, as the comparison between Baier's results for energy loss in a weakly interacting 
pion gas and in a weakly interacting quarki-gluon plasma would suggest. Certainly, the 
relationship between the energy loss of a hard parton traversing a medium and the energy 
density of the medium requires further elucidation and generalization. The perturbative 
expression for the radiative energy loss of an energetic parton, which is only known to 
leading order in the strong coupling a s , may get substantial corrections at higher order. 
A calculation of next-to-leading order corrections to the energy loss in perturbative QCD 
would be desirable. Although the energy of the penetrating parton provides a large scale, 
the momentum transfer of the scattering in the medium provides a second, much lower 
energy scale, which could enter into the NLO corrections. Also, the quantity q which 
provides the link between the observable AE and the energy density e is defined only in 
the context of the perturbative analysis. Since e itself is nonperturbatively well-defined, 
it would ultimately be desirable to understand the relation between parton energy loss 
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observables and e directly, possibly involving other nonperturbatively denned properties 
of the medium being probed by the hard parton. Keeping these caveats in mind, we 
nevertheless expect that the qualitative lessons encoded in Baier's results will survive in a 
more rigorous treatment. 

Jet quenching by itself can never provide a measure of v. It cannot differentiate 
between a weakly interacting quark-gluon plasma and a system whose gluon content is 
the same, but which has much less entropy density because the gluons are bound within 
hadrons and hence do not directly contribute to the entropy. In a hadron gas, the entropy is 
(roughly) a count of the hadrons, which are the thermodynamically independent degrees of 
freedom, whereas a hard parton "sees" the gluons within each hadron. A weakly interacting 
quark-gluon plasma and a hypothetical weakly interacting pion gas with the same energy 
density are equally effective at quenching jets, according to Baier's perturbative analysis, 
but the pion gas has a much lower entropy density and hence a lower v. Turning it around, 
if we imagine a hadron gas and a quark-gluon plasma with the same entropy density, the 
quark-gluon plasma has the smaller energy density, the smaller density of gluonic scatterers 
(only one per entropically active degree of freedom), and hence the smaller energy loss by 
a factor that scales like i/ 1 ^ 3 . The details of this calculation will differ depending on the 
precise nature of the composition of the medium, but the principle is general: Any hadronic 
medium will contain several gluons per hadron, but only the entire hadrons will contribute 
to the entropy. For a deconfined medium, on the other hand, each gluon contributes to the 
entropy individually, implying a smaller number of gluons (per unit of entropy) capable of 
scattering an energetic parton. 2 

Given a lower bound on s(to), in order to obtain a lower bound on v, and hence perhaps 
a demonstration of deconfinement, we will need to use jet quenching data to obtain an 
upper bound on e(to). Further theoretical work is required before this can be attempted 
quantitatively. Even if the relation between jet quenching data and energy density were 
in hand, current data may not be sufficient to provide an upper bound. The fact that 
the observed jet quenching in the most central collisions is quite large (the suppression 
ratio Raa *C 1) means that it may be difficult to derive an upper bound on q from these 
data [p7|| . The observed strong suppression implies that the observed hadrons are the 
leading partons from jets originating near the surface of the matter Q, as is also indicated 
by the absence (within current error bars) of any high-p^ hadrons from the away-side jet in 
the most central collisions pq| . Indeed, the currently ongoing careful study of the material 

2 As an aside, we should note that one proposed way of understanding the lattice QCD result that v in 
the quark-gluon plasma is somewhat less than that for an ideal quark-gluon plasma invokes the possibility 
that, in the energy density range relevant to the RHIC experiments, matter in the deconfined phase may 
contain colored bound states of quarks and gluons ]25| . This reduces the entropy density somewhat without 
reducing the density of gluonic scatterers in the medium. This is another illustration of our logic, but it is 
not directly relevant to the strategy for measuring v that we propose. Any bound on v extracted from data 
can be compared directly to lattice QCD calculations, regardless of the mechanism by which the strong 
interactions create the deviation of v from its ideal quark-gluon plasma value. Note also that colored bound 
state formation among quarks and gluons introduces a new contribution to the energy loss, that coming 
from ionization of the bound states . Neglecting this contribution in a future analysis of jet quenching 
would mean that the inferred energy density is greater than the true energy density, and thus would not 
interfere with the goal of obtaining an upper bound on the energy density. 
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recoiling against a jet originating from near the surface of the matter, for example via 
dihadron distributions, will provide us with data on new jet quenching observables and 
could teach us more about the matter produced in central collisions. However, since an 
upper bound on the energy density must take the form of a statement that if e were larger 
more suppression would have been seen, jet quenching data from noncentral collisions may 
be a more powerful source of information. Indeed, the away-side jet does not fully disappear 
in noncentral collisions, especially when the jets only have to fight their way through the 
hot matter in the narrow direction [28], giving hope that a quantitative upper limit for the 
value of q may not be far away. Data at higher px and with higher statistics is required, 
and should come from the Run-4 data set now being analyzed. 

Although the above paragraph is our true conclusion, we would be remiss to end 
without attempting to "plug in numbers just for fun", even absent a reliable upper bound 
on e(to). By applying the Bjorken argument to dEx/dy, both PHOBOS and PHENIX 
estimate that e(lfm/c)) > 5 GeV/fm~ 3 [El]. This lower bound on the energy density is of 
considerable interest in and of itself, even though by itself it cannot be used to constrain v. 
Indeed, if we take the lattice QCD calculation of v(T) as a given, the experimental lower 
bound on e tells us that at t = 1 fm/c the matter produced at RHIC is a quark-gluon 
plasma with a temperature T > 1.4 T c , well above the crossover. In other words, using 
(|l|), a lower bound on e, and the lattice calculation of u, we obtain a lower bound on T 
and a "demonstration" of deconfinement, albeit one that is unsatisfying because it uses 
the lattice calculation of v rather than testing it. In order to measure v and (presumably) 
demonstrate deconfinement, what is required is a lower bound on s(to) and an upper bound 
on e(to). If we adopt the conclusion from above that s(to) > 33 ± 3 fm -3 and, absent a 
reliable upper bound, suppose that e(lfm/c) is given by 5, 7 or 9 GeV/fm -3 , we would 
conclude that v > 71 ± 22, v > 26 ± 8 or v > 12 ± 4, respectively. A 5% determination 
of s(to) would reduce these error bars significantly, which motivates the theoretical effort 
needed to accomplish this goal. Another direction in which theoretical effort is needed is 
the modelling of the consequences of the variation of s and e across the transverse extent 
of the collision region, something we have not considered here. A stringent experimental 
investigation of low-p? event-by-event fluctuations is also required, in order to augment 
current theoretical evidence with experimental evidence against a strong first order phase 
transition, whose attendant entropy production would complicate the extraction of v that 
we propose. And, most important, these numbers make very clear the importance of further 
analysis of jet quenching theory and data with the goal of setting a reliable upper bound 
on e(io)- If there were experimental evidence that e(l fm/c) < 7GeV/fm -3 , this would be 
evidence for deconfinement. 

Looking further ahead, if RHIC data can provide interesting limits on the value of is, 
data from heavy ion collisions at the LHC should do even better. And, QCD predicts that 
if v is above the crossover at RHIC, its value will not increase significantly at the LHC. It 
is a greater challenge to devise a way of measuring v at lower energy densities using lower 
energy heavy ion collisions, where jet quenching is not observable. 
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